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Abstract 

The expansion of Kummer's hypergeometric function as a series of incomplete Gamma func- 
tions is discussed, for real values of the parameters and of the variable. The error performed 
approximating the Kummer function with a finite sum of Gammas is evaluated analytically. 
Bounds for it are derived, both pointwisely and uniformly in the variable; these character- 
ize the convergence rate of the series, both pointwisely and in appropriate sup norms. The 
same analysis shows that finite sums of very few Gammas are sufficiently close to the Kum- 
mer function. The combination of these results with the known approximation methods for 
the incomplete Gammas allows to construct upper and lower approximants for the Kummer 
function using only exponentials, real powers and rational functions. Illustrative examples are 
provided. 
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1 Introduction and preliminaries. 

Kummer's confluent hypergeometric function M(a, (3,x) = iFi{a, P,x) can be introduced in a 
number of equivalent ways [1] . For real values /3 > a > of the parameters and for all real a;, we 
can regard this function to be defined by 

M(a, (3, x) := — ^ r dt - tf'"-'e^' , (1.1) 

B[a,P-a) Jo 

where B is the Beta function (see the notes following this Introduction). The object of this paper is 
a series expansion for M derived from the above integral representation; for dealing more efficiently 
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with this expansion, we will work on the " reparametrized Kummer function" 



N{a,6,x):=[ dt t"''^ {1 - tf e'"'* for q e (0, +00) , 5 e (-1, +00), a; e [0, +00) . (1.2) 
Jo 

The knowledge of this function, with the variable x confined to the nonnegative semiaxis, allows to 
reconstruct the standard function M on the full real axis: in fact, for all a G (0, +00), /3 e (a, +00) 
and X G [0, +cxd) it is 

M{a,(3,-x) = — — ]- N(a,(3-a-l,x) and (1.3) 

B[a, p — a) 

M{a, f3, x) = —4 r e^Af(/3 -a,a-l,x) , 

B[a, p — a) 

where the first equality follows immediately from the definition of N, and the second one is inferred 

from the known identity [1] M {a, /3, x) = e^M{(3 — a, /?, —x). 

Let us fix the attention on N . If we expand in Taylor series the term (1 — tY in the integral 
representation (1.2) and integrate term by term, we obtain 

^(a.*,x) = g<^l(^, (1.4) 

where the 7's in the r.h.s. are incomplete Gamma functions, and ( );; is the Pochhammer symbol 
(see again the notes at the end of the Section; each summand in the above expansion also makes 
sense for a; = 0, for it admits a finite x 0"*" limit). A very simple case of Eq. (1.4) occurs when 
5 is integer: then (1 — ty is a polynomial of degree 5, and the above series for iV is a finite sum 
with fc = 0, 1, ...,5. Also, it should be noted that for integer a all the above 7's are elementary 
functions, expressible in terms of exponentials and powers of x. 

For large x, each incomplete function 7(0 + k,x) approaches the "complete" T(a + k); thus, one 
obtains a representation 

Nia,S,x)^"f^-^^-^ for. ^+00, (1.5) 

which has been extensively studied [5]; this must be dealt with carefully, because (for noninteger 
5) the series in the r.h.s. is not convergent, but only asymptotic in the Poincare sense [4]. 
In comparison with (1.5), the series (1.4) seems to have attracted less attention in the literature: for 
example, a recent survey [3] on computational methods for the Kummer function presents (1.4) as 
a "new method" for calculating this function. An argument of [3] proves convergence of the series 
for all x G [0, +oc); wc think that it is worth to continue the analysis of (1.4) deriving accurate 
estimates on the error for the series truncated at any finite order, and its rate of convergence. 
Obtaining these estimates is the aim of the present paper. Our approach will be based on some 
general inequalities about integrals of the form Jq dt ip(t)e~^*, of which Eq.(1.2) for N gives an 
example. The behaviour of these integrals for x — > +00 is the object of the classical Watson 
Lemma [4]; here we present a variant of Watson's analysis, allowing to derive upper and lower 
bounds on these integrals for all a; > by finite sums of incomplete Gammas. These bounds are 
almost self-evident, but in fact very useful in connection with N: after some technicalities, their 
application to N will yield the conclusions listed below. 

i) Upper and lower approximants of arbitrary accuracy can be derived for N, using finite sums of 
incomplete Gammas. 
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ii) The results i) imply convergence of the series (1.4) on the whole interval [0, +00), not only 
pointwisely but also uniformly in x; uniform convergence can be expressed using appropriate sup- 
norms. 

iii) The approximants for N obtained in i) are accurate on the whole semiaxis [0, +00) even when 
one sums very few terms (however, for very small x the Taylor expansion of N is generally more 
precise) . 

iv) The accuracy of the low order expansions in incomplete Gammas is essentially preserved if, 

for nonintcgcr a, one replaces these fimctions with some known approximants of Pade or other 
types. In conclusion, it is possible to obtain accurate upper and lower approximants for N for both 
integer and noninteger a using only exponentials, real powers of x and rational functions. 
Let us outline the organization of the paper. The rest of this Section fixes some notations and 
reviews some basic facts about Gamma functions and their approximations, the Beta function and 
the Pochhammer symbols employed extensively in the sequel. Sect. 2 presents the upper and lower 
bounds for a general integral of the Watson type in terms of incomplete Gammas. Sect. 3 analyses 
the error in the Taylor expansion of the funtion (1 — i)*, globally on the interval [0, 1); this is 
necessary for applying the bounds of Sect. 2 to the integral representation of N. Sect. 4 contains the 
main results about the expansion of N by incomplete Gammas, and the error for the expansion 
truncated at any finite order. Bounds for the error are derived analytically, both pointwisely and 
uniformly in x; some numerical tests on these theoretical bounds are presented, for an appreciation 
of their reliability. Sect. 5 illustrates by some examples the low order expansions for N ^ and the 
elementary substitutes derived from them replacing the incomplete Gammas with their known 
approximants by exponentials, powers of x and rational functions. 

Miscellaneous notations. Throughout the paper N and Ng are, respectively, the nonnegative 
and the positive integers. The symbol ^ indicates that two functions or sequences are asymptotic 
in the elementary sense, i.e., that their ratio tends to unity. 

On Gamma functions. We use the standard notations for the incomplete and "complete" 
Gamma functions; in the ranges indicated below, we can regard them to be defined by 

px r+00 

^{v,x) := / ds s^-^e-^, T{v) := / ds s^-'^e'^ for v e (0, +00), a; e [0,-|-oo). (1.6) 
Jo Jo 

For either v ox x fixed, the incomplete function has the asymptotics 

, X f ~ x" Iv for X 0"*", „, 



- x'^/v 


for 


X - 


-0+, 




for 


X - 


+00 ; 




for 








V - 


+00 . 



Of course. 



(1.8) 



j^^{u,x)=x''-^e-- ; (1.9) 



also, we have the recursion rule 

7(1/, x) = {u- 1)7(1/ -l,x)- x^-^e-"" . (1.10) 

Concerning the complete function, we will continuously use the relations 

V{v+\) = vV{v) (1.11) 

and r(n+ 1) = n! for n e N. Sometimes we need the standard extension of the complete function 
r from [0, +00) to R\ (— N), which is uniquely determined asking (1.11) to hold everywhere on 
this enlarged domain. It is V{y) ^ 00 for z/ — > p e — N. 
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Practical evaluation of 7(1^, x) for fixed v G (0, +cx)). This is very simple for integer v: 
using the recursion rule (1.10), we reduce the computation to the case 

7(l,ar) = l-e-- . (1.12) 

For arbitrary u one constructs approximants of '){y,x) of different types, for x small or large. 

Some elementary approximants are derived from Taylor's formula applied to the exponential e~* 
in Eq.(1.6). For s > and m G Nq, we have = X)J1^^(— l)-'s-'7i! + (— l)'"e~'^s'"/m!, where 
c e [0, s]; of course, < e~° < 1 if s € [0,x]. Thus 

V t}ls^ + < e- < y t}^s^ + ElM.m (113) 

1! m! V. ml 

j=o •' j=o •' 

for X e [0, +00), s e [0,a;], m e Nq, where 

..(.) := min ((-1)™, (-l)-e-^) = { ^ ^^^j' (1.14) 

- -ax {i-ir, i-l)-e-^) . { - (1.15) 

We substitute the inequalities (1.13) in the definition (1.6) for 7, and integrate term by term. This 
gives 

Tm(z/, x) < 7(1/, x) < Tto(i/, x) for u e (0, +00), X e [0, +00), m e No , (1-16) 

(-l)J x'^+^ , ^„(x) 
T„(i/,a;):=> — — + ■ ; , (1.17) 

The "Taylor" approximants r™, are interesting for small x only; for large a;, a completely 

different approach must be employed. As a preliminary, one observes that repeated application 
of the recursion rule (1.10) starting from a noninteger value of v reduces the problem to the case 
V e (0, 1): in this case reliable approximants are known for large x, mainly of two types. 
Firstly, there is a sequence of "Laurent" approximants, involving negative powers of x. These 
are constructed working on the representation 7(1^,3;) = T{y) — ds s''~^e~*; by repeated 
integrations by parts, one infers [4] 

T{v) - x^'-^e-'^kqiv, x) < 7(1/, x) < r(i/) - x^-^e-'^Xgiiy, x) (1.18) 

for all 1/ e (0, 1), X e (0, +00) and q £ Nq, where 

fe=0 k=0 

(m)o:=1, (m)/. := Mm -!)(/« -2)...(m-/i+ 1) for /i S Nq, (1.20) 
X+ := max(x, 0) , x~ := min(x, 0) for x S R. 
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Typically, these Laurent approximants are accurate for very large x but less satisfactory for in- 
termediate values, close to 1. This defect is overcome by the Pade approximants, constructed by 
more refined techniques [2] : these have the form 

T{v) - x''-^e-'=\lq{v, x) < 7(!/, x) < T{i^) - x''-^e-'''Kq{v, x) (1.21) 

for ^ e (0, 1), a; € (0, +oo) and g S N, where both Ilg, TTg are ratios of polynomials of degree q in 
x; the lowest order cases are 

X X -\- \ 

'Ko{v,x) :=0, Iio{v,x) := 1 ; ■ki{v,x) := — — , Ili{u,x) := — — ; (1.22) 

X + 1 — v X + 2 — u 

, , x{x + 3-i^) TT / ^ + (5 - z/).T + 2 

Tr2{iy,x) := ^ — — — -, U2{u,x)~ 



a;2 + 2(2 - iy)x + (1 - p){2 - v)' ' ' ' a;2 + 2(3 - v)x + (2 - z/)(3 - v) ' 

Of course, one can match the approximants of the above different types to get globally accurate 
estimates. For example, let us evaluate 7(1/2, x) using the Taylor bounds of order m = 4, and the 
Pade bounds of order g = 1; this gives the inequalities h{x) < 7(1/2, x) < H{x) where h is the 
maximum of the two lower bounds, and H the minimum of the two upper bounds {h is the Taylor 
lower bound for < .t < 1.16, and the Pade one for x > 1.17; H is the Taylor upper bound for 
< a; < 1.48, and the Pade one for x > 1.49). The relative uncertainty {H-h)/{H + h) is < 0.005 
for all X e [0, +00), and attains its maximum at a point xq € (1.48, 1.49). 
Pochhammer's symbol. For a e R and fc e N, this is defined setting 

(a)o:=l, (a)fe :=a(a + l)(a + 2)...(a + fc-l) ; (1.23) 

one has 

(a)fe+i = {aUa + k) = a{a + 1). , ^^±i^ = E ^ (1-24) 

fe=0 

(the first two equalities are elementary, the last one is proved recursively). It must be noted that 



F Q: + k) 

Wk = X , (1-25) 

r(a) 

(provided that, for a € — N, the r.h.s. be intended as a limit from nonintegcr values; incidentally, 
the above relation allows to compute r(l/2 + k) starting from F(l/2) = ^/tt). The behaviour of 
Pochhammer's symbol for real a and large fc e N is described by the equation 

, , / ~ A;! fc"-Vr(a) for a ^ -N, fc ^ +00 , 

(")'=\=0 forae-N,fc>l-a; ^^"^^^ 

(the statement for a G — N is evident, the other one follows from known properties of the Gamma 

function [4]). 

The Beta function B. The standard definition reads [1] 



B{n,v):=[ dtt^'-'^il-t)''-'^ for e (0,+oo); 

Jo 



(1.27) 



this function is related to the complete Gamma by the identity 
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Elementary facts about N. Throughout the paper we stick to the definition (1.2) of N{a, 5, x), 
for all a € (0,+oo), S G (— l,+oo) and x G [0,+oo). Even though our main interest is the 
representation (1.4) for N in terms of incomplete Gammas, it is convenient to review the Taylor 
expansion of N that will be sometimes compared with (1.4). The derivation of the Taylor expansion 
is elementary: we substitute in Eq. (1.2) the inequalities (1.13), with s = tx, and we integrate 
term by term. Recalling the definition (1.27) of B, we obtain the final result 

tm{a, S, x) < N{a, 6, x) < T^{a, 5, x) (1.29) 

for a e (0, +oo), 5 e (—1, +oo), x e [0, +oo) and m € Nq, where 

tm{a, 6, x) := V ^^^fi(a + k,d+ l)x'' + ^^^^S(a + m,6+ l)x"' , (1.30) 

r„(a, S, x) := y ^^B{a + k,S+ l)x'= + ^^^B{a + m,S+ 1)^™, 
f-^ k\ ml 

k=0 

and Um, Um are as in Eq.s (1.14-1.15). Sending m to +oo in Eq. (1.29), we obtain the power series 
representation 

N{a,S,x) =y^^B{a + k,S+l)x'' , (1.31) 

fe=0 

converging for all x in the chosen domain [0, +oo). The power series expansion is often used as a 
definition of the Kummcr function, alternative to the integral representation adopted here as our 
starting point. As a final remark about (1.31), we note that Eq.s (1.28) (1.25) imply 

B{a + k,S+l)= B{a,S + l) . (1.32) 

[a + d + 1 fe 



2 Integrals of the Watson type. 

By a Watson type integral, we mean an (improper) Riemann integral depending on a real parameter 
X, of the form 

W{x):= [ dt i;{t) e-''* , b e {0,+oo], tj; G C{{0,b),R) ; (2.1) 

in the sequel, the expression "W^(x) exists" will mean that the integral (2.1) is convergent, for 
the considered value of the parameter x. Watson's Lemma [4] was originally formulated as a 
statement on the behaviour of W{x) for x +oo, under appropriate conditions on ^. Here we 
present a simple Watson- type statement giving upper or lower bounds on W{x), without sending 
X to inRnity: this is an inequality, holding for both small and large x. 

2.1 Proposition. Assume 

n 

^{t) > {resp. <)J2Pkt^'' e (0, 6), (2.2) 

fc=0 

where n G N, pk G R, Xh & (—1, +oo) for k = 1, n. If x G (0, +oo) and W{x) exists, then 

Wix) > (resp. <) Ep.%tilM . (2.3) 

fe=0 
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In the above equation, for b = +00 one intends bx := +00, 7(Afe + 1, +00) := r(Afe + 1). 

Ifb < +00 and W(0) exists, the inequality (2.3) holds for it replacing each term 7(Afe + l, bx)/x^''^^ 

with its 3; ^ 0+ limit, which equals h^''^^ / {\k + 1). 

Proof. Wc simply insert the inequality (2.2) in the definition (2.1) of W{x); this gives W{x) > 
(resp. <) Y^^=oP^ lo t^''e~^^. For x € (0, +00) and b either finite or infinite, the variable change 
t = s/x in the integrals yields the thesis. 

For a; = 0, and h finite, each integral in the previous sum equals b'^''~^^/{Xk + 1), which is the 
a; ^ 0+ limit of 7(Afc + 1, bx)/x^''+^ due to (1.7). o 



3 The error in Taylor's expansion of (1 — t)^. 

The function N(a,(3,x) of Eq.(1.2) is a Watson type integral of the form (2.1), with 6=1 and 

V'(t) = V(a,<5,t) :=t"-i(l -t)*. 

We want to find both upper and lower bounds for in terms of powers of t, of the form (2.2); the 
problem can be reduced to studying the Taylor expansion about t = of the function 

t e [0, 1) (1 - i)* e C°°([o,i),R) . (3.1) 

For our purposes, it is essential to discuss carefully the behaviour of the error everywhere on the 
interval t e [0, 1); this is the object of the following 

3.1 Lemma. For all S G R and n G Nq, there is a unique function Pn{S, •) € C°°([0, 1), R) such 
that 

It is 

p^{d, t) = / du {1- u)""^ (1 - to)''-" , sign pn{6, t) = const. = sign (-(5)„ . (3.3) 

(n-1)! Jq 



(l-tf = J2^^t'' + pn{6,t)t'^ V<G[0,1). (3.2) 

fe=0 



Furthermore, 



-(l-<5)„_i/(n-l)! if S>0, 
Pn{S,0) = ; .lim p„(J,i) ={ if S = 0, (3.4) 

f oo if 6 <0. 



For all t € [0, 1), the derivative of pn{S,») is described by 

p'niS, t) = /' du u{l - n)"-i (1 - tu)'—' , (3.5) 

sign p'j^{6,t) = const. = sign (— ^)„+i . 
(All statements on the signs include the zero case, intending sign{0) := 0.^ 

Proof. Taylor's formula for an arbitrary / G C""([0, 6), R), (m > n) ensures the existence of a 
unique function p„ G C"'-"([0, b), R) such that 

n — 1 

= E fc! /^'H0)t^ + /)„(i)r for t e [0, b) , (3.6) 

fe=0 
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which has the integral representation 



The function /(t) = f{5, t) = {l- t)^ on [0, 1) has derivatives 

f^'\t) = ^ (1 - = - ^)'"' ; (3-8) 

this yields Eq.s (3.2), (3.3). The integral in Eq.(3.3) is clearly positive, so pn{S,t) has the sign of 
the coefficient {—S)n- 

In general, for the function /3„ of the expansion (3.6) we have Pn(0) = /^"■'(0)/n!, which equals 
(— ^)ra/^! in the present case. The results (3.4) for lim^^i- pn{S,t) can be derived from (3.2): 
if 5 > 0, the t 1~ limit computed from (3.2) is in fact — X]fc=o (^^)fc/''''' but this equals 
— (1 — S)n-i/{n — 1)1 due to the third relation (1.24). (This case could be treated alternatively by 
the integral representation (3.3), which implies 

{-S)n f' _ {-S)n 1 _ (1 - S)n-1 



lim Pni5,t) = - — [ du {1- u) 

t^i- [n - 1)! Jo 



t^i-^"^'^ („_1)! ^ ' {n-iy.S (n-1)! 

due to the second relation (1.24)). 

To obtain the representation (3.5) for p'^{5,t), one dcrivates Eq.(3.3) and then employs the first 
identity (1.24) to express the coefficient before the integral in the form {—6)n+i/{n — 1)!. The 
integral in (3.5) is positive, so the derivative p'^ has the sign of {—5)n+i- o 

Remark. Pn{S,») and its derivative are hypergeometric functions of the Gaussian type [1]. The 
Gaussian hypergeometric function F = 2F1 with parameters a, b, c can be defined setting 



/ du u^-\l - uY-''-\l 

,c-b) Jq 



F{a,b,c:t) -.^ — - I duu''-\l-uy-''-'il-tu)-'' . (3.9) 
By comparison, wc see that 

p45,t) = ^^ F{n-S,l,n+l;t) , p^S, t) = i^(n + 1 - ^, 2, n + 2; t) . (3.10) 

The previous Lemma implies 

3.2 Corollary. For 6 & R and n G Nq, consider again the function pn{S,») of Eq.(3.2); then, 
its sup and inf are as follows. 
For 5>Q, 

inf pr,{5,t) =rn{5), sup Pn{5,t) = Rn{S) , (3.11) 

*6[o,i) te[o,i) 

r ix\ ._ f (-^)» _O^^^S)n-i\ _ / the left term if (-5)„+i > 0, , , 

V n\ ' (n-1)! J - \ the right term if {-6)r,+i < 0, ^'^■'^> 

^' \ n\ ' (n-1)! J {the left term if{-S)n+i<0- ^ ' 
For 6 < 0, it is 

inf pn{S,t) = rn{S) , sup pn{d,t) = +00 , (3-14) 

*e[o,i) t6[o,i) 



with 



r„(5) := . (3.15) 
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Proof. In any case pn{S,») is monotonic, so its sup and inf are the limit values at f = or 
t = 1. The thesis follows from the results of the previous Lemma about these limits and the sign 

oip'niS,')- O 

Remarks, i) For all n E Nq, it is 

-^^^ = 1^^ if^eR\{0}; -(i^%i = -l, t^ = if^ = (3.16) 

(n— 1)! on! (n - 1)! n! 

(the first relation above follows from the second identity (1.24)). 

ii) Many of the previous statements involve the signs of the Pochhammer symbols of —S; therefore, 
the following tables can be useful. For a real noninteger S, with integer part [S\ , and for all e N 
it is 

r (-1)'=^ if S>0,0<k< [5], 

sign(-5)fc=<^ (-1)1^1+1 if S>0,k>[6] + 1, (3.17) 
[ 1 if 5<0. 

For S a relative integer and fc e N, 

if 6>0,0<k<6, 
sign(-5)fe = { if (5 > 0, > (5 + 1, (3.18) 

if 6<0. 

(Of course, the vanishing of (— ^)fe for ^ e N and fc > ^ + 1 reflects the fact that, in this case, 
(1 — ty is a polynomial of degree S). o 




4 Expanding the Kummer function. 

We are now ready to derive the main results on the above expansion. As in Sect.l, we fix the at- 
tention on the reparametrized function N{a, S, x) of Eq.(1.2) for a G (0, +oo), S G (—1, +oo), x G 
[0, +00) (the relations of N with the standard Kummer function M have been exploited in the 
Introduction). 

We first discuss the case (5 > 0, where the function pn{S, •) of the previous Section is bounded both 

from below and above for all n. 

4.1 Proposition. Let 6 G [0, +00), a G (0, +00) and n G Nq. For all x G [0, +00), it is 

gn{a, 5, x) < N{a, 6, x) < Gn{a, S, x) , (4.1) 

where 

..(a.*..):=|:H)i2(£±^ + ,„(,)2(;i±^:£), (4,2) 

fe=0 

tr. X ^\ V" {-S)k"/{a + k,x) 'y{a + n,x) 

Gn{a,5,x) ■■=2^—^ — +Rn{5) 

fc=0 

and rn, Rn are as in (3. 12-3. IS). For x = 0, each term 7(a + h,x)/x°'~^'^ in the above is intended 
to mean its a; ^ 0+ limit, equal to l/{a + h). 
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Proof. As already noted, N{a,S,x) is a Watson type integral (2.1), with 6=1 and = 
ip{a,6;t) := t"~^{l — f)*. On the other hand, Corollary 3.2 gives the inequalities 



fe=0 fe=0 



whence 



|j Likf«-i+fe + < _ < J2 L^e-^+f + R^{S) ; (4.4) 

fc=0 ■ k=0 

thus, the thesis follows immediately from Prop.2.1. o 

We pass to the case (5 < 0, where the function p„((5, •) of the previous Section is bounded only 
from below; in this case, a different method will be used to get an upper bound on N{a, 6, x). 

4.2 Proposition. Let 6 G (—1,0), a G (0, +oo) and n G Nq. For all x G [0,+oo), it is 

gn{a, S, x) < N{a, S, x) < Gn{a, S, x) . (4.5) 
Here g„ is as in Eq. (4-2), with r„(5) = (— 5)„/n! as prescribed by (3.15); furthermore 

Gn{a,6,x):=}^— ^^^^ +Sn{a,5) ^^^^ , (4.6) 

fe=0 

S^ia, S) := ^ J . (4.7) 

(As usually, for x = one intends the previous definitions in a limit sense). 
Proof. In this case Corollary 3.2 implies 

E + r„((5) < t«-i(l - tf . (4.8) 

fe=0 

Using again Prop.2.1 on Watson integrals wc infer from here the lower bound g„, of the form (4.2). 
In order to obtain the upper bound, wc use a known recurrence relation for the Kummer function. 
This is the identity M(a, /3, x) = M{a, x) + {ax/ P{P+l))M{a+l, /3 + 2, x) [1] that becomes, 
in terms of the reparametrized function N, 

N{a, S, x) = , T N{a, 6+l,x)- ^— 7V(a +l,S+l,x) . (4.9) 
d + 1 + 1 

The identity holds for all S E (— l,+oo) and a. E (0,+oo); with the assumption S E (—1,0), the 
N functions in the r.h.s. depend on the parameter 6 + 1 G (0,1) and their expansions in terms of 
Gamma functions can be performed via Prop.4.1. More precisely, we have 

N{a,S + l,x) <^ __^+iJ„(5 + i)____ , 

, <- , N v-^ (— (5 — l)fc 7(a + fc + 1, x) ,^ ,,7(a + n + l,x) 
N{a + l,5+l,x)>Y^ fc, +rn{5 + , (4.IO) 

fe=0 
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Rn[0 + i) = j , r„(d + l)- 



(n-1)!- 

Inserting these inequalities into the identity (4.9), we get an upper bound for N{a,d,x). This 
upper bound can be written as in Eq.s (4.6-4.7) reexpressing j{a + k + l,x) and 7(a + n + l,x) 
by means of the recursion rule (1.10), and then manipulating some of the occurring Pochhammer's 
symbols via the identities (1.24). <> 

Remarks, i) For 6 e (-1,0) and aU a € (0,+oo), it is < r„(5) < Sn{a,6). 

ii) The remaining part of the Section will be devoted to the n — > +oo limit in Prop.s 4.1, 4.2, 
yielding the convergent series expansion (1.4) for N. 

iii) Comparing the definitions (1.2), (1.27) for A'' and B, we see that 

N{a,S,0) = B{a,S + l) . (4.11) 

Thus, Prop.s 4.1, 4.2 imply 

^k\a + ka + n~^ ' ~ f-^ k\ a + k a + n ^ ' 

fe=0 fe=0 

For n — > +00, this gives a series expansion for the Beta function. o 
The error €„ in the expansion of JV. For conveniency, let us define 

e„(a, <5, X) := iV(a, ^, x) g ^'fc^' "'^'^^H-t' ^ (4-13) 

/c=0 

this is the error in the approximation of A'^ by means of n incomplete Gammas. Prop.s 4.1 and 4.2 
can be rephrased as 

e„(a, (5, x) < e„(a, ^, x) < En{a, S, x) (4-14) 
for a e (0, +oo), 6 G (—1, +oo), x G [0, +oo), n G No, where 

en{a,d,x) := rn{d)^^^^^^^ , E^ia,d,x) := (i?„(^) or Sn{a,5)) ^^'^^t+n''^ • (4.15) 

The coefficient in the definition of En is Rn{^) for S E [0, +oo), and Sn{a,S) for S G (—1,0); as 
usually, 7(1/, x)/x'^ := l/v for a; = 0. The quantities e„. En are the theoretical bounds on the error 
e„, derived from the previous Propositions. Of course, by the asymptotics (1.7) of 7, fixing n we 
have 

en{a,S,x) ~ rn{S)^^^^^^^,En{a,S,x) ~ (i?„(5) or Sn{a,S)) ^^^^ for x ^ +00 (4.16) 

(to get a better insight into e„ and E^ as functions of x, one could employ for 7(0 + n,x) the 
approximants described in Sect.l, for example the Taylor approximants for small x and the Pade 

ones for large x). 

Let us discuss the error bounds when n is large, starting from a pointwise analysis in x (a uniform 
analysis is be performed in the sequel). 

4.3 Proposition. Fix a G (0, +00), 6 as below and x G [0, +00). Then 

er-^/Sr{-S)n^-^^ for S € (-1, +00) \ N, [S] even, n +00, 
en{a, S,x){ - e-^/r(-5)n2+'5 for S e {-1, +00) \ N, [6] odd, n +00, (4.17) 
forSGN,n>S+l; 
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En{a,S,x) < 



' ~ e'^ lV{~S)r?^^ for b G (0, +ocO \ N, \S\ even, n +00, 

~ e-''/ST{-6)n^+'^ for S G (0, +00) \ N, [6] odd, n +00, , . 

~ e-^/(l + (5)r(-(5)ni+'' /or (5 G (-1, 0), n ^ +00, ^ ' 

= for 5 e N, n>5+l. 



In any case e„(a,5, a;), En{a,5,x) — > 0, and thus 

fc=0 

(in the special case S € N, the sum in the l.h.s. equals N{a,l3,x) if n — I > 6). 

Proof. The definitions of the coefEcients r„, -R„ and Sn, with Eq.s (3.17-3.18) and (1.26) on the 
signs and asymptotics of the Pochhammer symbols, yield the conclusions 

-l/r(l - 6)n^ = l/Sr{-6)n^ S G (-1, +00) \ N, [d] even, n +00, 
rn{S){ ~ l/r(-^)ni+'5 (5g (-l,+oo)\N, [5] odd, n^+00, (4.20) 

(5eN, n>(5+l; 

- l/r{-S)n^+^ 6 G (0, +00) \ N, [6] even, n +00, 

Rn{S){ ~ -l/r(l- J)n'' = l/5r(-(5)n'' 5 G (0, +00) \ N, [5] odd, n ^ +00, (4.21) 
5gN, n>(5+l; 

(5TT)fRK' '^e(-l,0),n-+oo. (4.22) 

Furthermore, Eq.(1.8) implies 7(0: + n, x)/x"~^" e~^/n for n +00; the thesis follows immedi- 
ately, o 

Of course, the absolute value of the error e„ can be bounded as well by means of Prop.s 4.1, 4.2, 
which imply 

\e^{a,6,x)\<£n{a,S,x) := X„(a, (5)I^^±^ , (4.23) 

for all a G (0, +00), 6 G (—1, +00), x G [0, +00), where 

/ max(|(-<5)„|/n!,|(l-(5)„_i|/(n-l)!) for 5 G [0, +cx3) , 
' ■ [ S'„(a, i5) for5G(— 1,0). ^ ' ' 

The analogue of Prop. 4. 3 for the theoretical bound is the following 

4.4 Proposition. For fixed a G (0, -hoo), 5 as below and x G [0, -hoo), it is 

( ~ e-^/|(5r(-(5)|ni+^ for S € {0, +00) \ N, n ^ +00, 

£n{a,6,x)i ~ e-^/(l -h 5)r(-(5)ni+'' /or (5 G (-1, 0), n ^ -hoc, (4.25) 
[ = for5GN,n>5 + l. 

Proof. The thesis follows from the known behaviour of j{a + n, x)/x"+" and from the asymptotics 

- l/\5r{-S)\n^ for 6 G (0, +00) \ N, n ^ -|-oo, 

Xn{a,6){ l/{l + d)T{-6)n^ for (5 G (-1, 0), n ^ +00, (4.26) 
= for 5 G N, n > (5-M. 
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Reliability of the error bound |en| < £n' numerical tests. Of course, one would like the 

actual error |e„| to be close to its theoretical estimator £„■ Figures 1-3 report the ratio |e„|/5n for 
fixed values of a, d,x as a function of n. They were generated by the MATHEMATICA package, 
using its internal routines for computing numerically the Kummer function and the incomplete 
Gammas. For reasons related to numerical accuracy, the analysis has been limited to 1 < n < 13; 
the outcomes seem to indicate that |e„|/f„ approaches for large n a constant value not far from 1, 
which was just the hoped result. For a better appreciation of the trends, the points in the pictures 
corresponding to n = 1, 2, 13 have been interpolated by continuous or dashed lines, according 
to the considered values of x. 
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" 1 3 5 7 9 11 13 

Figure 1: Ratios |e„|/f„ for a = 2, (5 = 3/2, x = 1 (continuous line) and x = 6 (dashed line). 
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1 3 5 7 9 11 13 

Figure 2: Ratios \en\/£n for a = 3, (5 = 15/2, x = 1 (continuous line) and x = 6 (dashed line). 
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0.4 
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" 1 3 5 7 9 11 13 

Figure 3: Ratios \en\/£n for a = 2, S = —1/2, x = 1 (continuous line) and a; = 6 (dashed line). 
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Uniform analysis of the error. We consider again the error e„ defined by (4.13), and analyse 
it uniformly in x] to this purpose, it is convenient to introduce an appropriate functional setting. 

For (7 G [0, +oo), we consider the space 

:={/gC([0,+oo),R) I +oo> sup x''\f{x)\:=\\f\\,}; (4.27) 

a:e(0,+c!o) 

this is a Banach space with || ||(^ as a norm. We note that 

^^^^'^ e C for > 0, a e [0, v] , iV(a, (5, •) G for a > 0, 5 > -1, a G [0, a] (4.28) 

(in the case of 7(1/, •)/»'^, this follows from the asymptotics (1.7); in the case of N{a,5,») the 
statement follows from Prop.s 4.1, 4.2 with n = 1, and from the asymptotics of the Gammas 
therein) . 

The following result will be used to bind the cr-th norm of 7(1^, •)/•" ■ 

4.5 Lemma. For all a,x G [0, +00) and fi G (0, +00), it is 

< :i^^^±^ < ^ (4.29) 
~ xi^ ~ IJ. ^ ^ 

(intending this in a limit sense for vanishing a or x: in particular, a" := 1 if a = 0). 

Proof. The lower bound is obvious, let us prove the upper bound. We can assume a > 0, and then 
recover the zero case taking the limit. So, let tr > 0; due to (1.7), the function x 7(/x + a, x)/x^ 
vanishes for x approaching both 0+ and +00, and thus it attains its absolute maximum at a point 
c G (0, +00); c is a stationary point, i.e., a solution of the equation 

= 4- 2^t±^ = i f x-e- - ^lil^l^) . (4.30) 

dx xt^ X \ xf^ J ^ ' 

Thus, the maximum of the function is 

7(/^ + tT,c) ^ £-£2 < ^"e-" , 

ci^ M ~ /i ' ^ ' ' 

the first equality follows from Eq.(4.30) and the subsequent inequality is elementary, since the 
function x G [0, +00) x^e~^ attains its maximum at a; = u. o 

Remark. There is numerical evidence that the maximum point c = c^^- of the above function is 
unique, and that 'yi^jJi + (J, c^a) / {c^aY is not only bounded by a'^e~'^ /fi, but in fact asymptotic to 
the latter for ji — » +00. <> 

4.6 Corollary. For all a G [0, +00) and v G (cr, +00), it is 



< • (4.32) 

ly — (7 

Proof. The above norm equals sup3,£[o,+oo) l{'^^x)/x'^~'^; now apply the previous Lemma. o 
4.7 Proposition. Let a G (0, +00), 5 G (—1, +00), a G [0, a]. For all n G iVo it is 



||e„(a,^,-)L < Xr,{a,5) " ^ , (4.33) 
a — (7 + n 



with Xn is as in Eq.(4.24)- 
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Proof. Eq.(4.23) implies 



7(a + n, •) 



\K{a,6,*)\\^<X^{a,S) ^^^^ 
now, use the previous Corollary. 

4.8 Proposition. Let a € (0, +00), 6 £ (— 1,+cxd), a e [0,a]. Then 



(4.34) 



r.h.s. of (4.33) 




1+5 



'/m-s)\n 





a+<5 



for S G (0, +00) \ N and n — » +00 
for 6 € (—1,0) and n +00 
for S e N and n>S+l; 



thus 



E 

fe=0 



i-d)k 7(a + fc,») 



N{a, 6, •) for n — > +00, in the Banach space C„. 



(4.35) 
(4.36) 



Proof. Use the asymptotics (4.26) for X„. 



The accuracy of the expansion of N via n incomplete Gammas can be discussed from a slightly 
different viewpoint, where the attention passes from the "absolute error" e„ to the relative uncer- 
tainty 

Gn{,a,5,x) - gn{a,5,x) 



^n{oi-,5,x) := 



(4.37) 



\Gn{a,5,x) +gn{a,5, x)\ 

here G„ and gn denote, as usually, the upper and lower bounds for TV in Prop.s 4.1, 4.2 (^). From 

the expressions of these bounds, it is clear that ^n(a, 8, x) = 0(1/ x"") for x — > +00. 

In the next section, the a;-dependence of ^„ will be appreciated numerically in a number of examples. 



5 Low order approximants for AT. 

Let us show by some examples that Prop.s 4.1, 4.2 give good approximants for N even if they arc 
used with low n. The second example will show that a good accuracy is kept even if we replace the 
incomplete Gammas with the approximants in terms of exponentials, real powers of x and rational 
functions reviewed in Sect.l (in the other examples the Gammas are elementary functions, with 
exact expressions in terms of exponentials and powers of x). In all the examples, the low order 
approximants for A'' will be compared with the Taylor ones (also presented in Sect.l), which are 
found to be more precise for small x, namely, for x ranging from to some value of the order of 
unity. 

i) We consider the case a = 2, 6 = 3/2 and apply Prop. 4.1 with n = 2. This gives 

g{x) < N{2, ^, x) < G{x) for all x G [0, +00) , (5.1) 

where 

n(^^ - ^(2^^) 37(3, x) ^ 37(4, x) _ 7(2, x) 37(3, x) 7(4, x) 



1 A comment is worthy on the absolute value in the denominator of ^n- The function N is positive for all a, S 
and X 6 [0, +00), so its upper bound G„ is also positive. The lower bound gn and the sum can be negative 

in some peculiar cases with large a, 6 and very small n, x. For example, Gi (10, 10, 0) + gi (10, 10, 0) < 0. 
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9(0) G(0).l 

All the above incomplete Gammas are of integer order, so they can be computed by the recursion 
method recalled in Sect.l: 7(2,2;) = 1 - {x + l)e-'=, 7(8,2;) = 2 - (a;^ + 2a; + 2)e~"', 7(4, a;) = 
6 — (a;^ + 3a;^ + 6x + 6)e~^. The conclusion is 

. X _ 1 3 9 / 1 7 3 9 \ _^ 
^^'''-x^ a;3+4a;4 + Ua; + 8a;2+4a;3 4a;4j'' ' ^^■'^> 



^/ N 1 3 3 



(—-- 
V2a;2 a;4 



Figure 4 contains the graphs of g, G defined as above, and of N{2, 3/2, •) as computed numerically 
by the MATHEMATICA package; the chosen interval for x is [0, 7]. For a better appreciation of 
the separation between N{2, 3/2, •) and the bounds g, G, Figures 5 and 6 report the graphs of the 

same functions for x G [0, 1] and x G [6, 7], respectively. 

Let us describe the discrepancy of the bounds g, G through the relative uncertainty (^) 

= Ira ■ 

This function is decreasing: ^(0) = 1/7 < 0.15, ^(1) < 0.098, ^(3) < 0.045, ^(7) < 0.012 and 
^(a;) — > for a; — > +00. 

The approximants for A^(2,3/2, a;) constructed above can be improved for small x, using in place 
of g, G the Taylor bounds (1.29-1.30); with the choice to = 4, these are 

t{x) < N{2, ^, x) < T{x) for x e [0, +00), (5.5) 

4 16 16 2 128 3 64 4 _^ 

t(x) := X H x^ x'^ H x^e ^ , (5.6) 

^ ^ 35 315 1155 45045 135135 ' ^ ' 

^. , 4 16 16 2 128 3 64 4 

T(x) := X H x^ x-' H x^ . 

^ ' 35 315 1155 45045 135135 

The final bounds arising from (5.1) and (5.5) are 

j{x) :=max{t{x),g{x)) < N{2,^,x) < J{x) := mm{T{x),G{x)) (5.7) 

for all X e [0, +00). One finds numerically that: j{x) = t{x) for < x < 1.92 and j{x) = g{x) 
for X > 1.93; J{x) = T{x) for < a; < 2.16 and J{x) = G{x) for x > 2.17. Let us introduce the 
relative uncertainty 

Then: ?7(0) — 0, 77(1) < 0.0021; sup^^[Q^^^-^r]{x) < 0.062, the sup being attained at the point 
a;o G (2.16, 2.17) where T{xo) = G{xo): ri{x) = £,{x) for x > 2.17. o 

ii) We keep d = 3/2 as before, but assume a = 1/2. Prop. 4.1 with n = 2 gives 

g{x) < N{^,^,x) <G{x) for all a; e [0, +00) , (5.9) 



^here and in the forthcoming examples, the bounds for which we evaluate the relative uncertainty are always 
positive, for all x S [0, +00). 
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where 

n(^^.-l^dllA 37(3/2, a;) 37(5/2, a;) . ._ tOM 37(3/2, 7(5/2, x) 

^^"'^ - 2^^ + 8a;5/2 ' - ^^37^ + 2a;5/2 ^^"^"^ 

23 fi 
.(0) = ^, G(0)^- 

One can express 7(5/2, •), and 7(3/2, •) in terms of 7(1/2, •) by means of the recursion rule (1.10); 
in this way, the previous bounds take the form 

^/x/l 3 3\,1,/1 3 

From the above result we can infer upper and lower bounds for 7V(l/2, 3/2, x) in terms of elementary 
functions, using for 7(1/2,2;) the Pade approximants described in Sect.l. 

Let us recall the general Pade approximants (1.21-1.22). To get a lower bound for g{x), we 
substitute in Eq.(5.11) the inequalities (l/Vi + 9/32 :rr^/'^yil/2,x) > (l/Vi+9/32 x-^^'^)x 
(the Pade lower bound of order q = 1), —3/4 a;^'^/^7(l/2, x) > —3/4 x~'^/-^x (the Pade upper 
bound of order 1); we proceed similarly to obtain an upper bound for G{x), using again the 
approximants with q = 1. In this way, for all x G (0, +00) we find 

9ix)>pix) + ^ + + 



3 2 9 9 9 



2x(l + 2x) .T(3 + 2a;) 16a;2 16a;2(3 + 2x) 16a;3(3 + 2a;). 

GixXPix) .= ^-i^ + i^ + f 1-^+ 

L.[x) s n^) . ^ ^^3/2 + g^5/2 + V a; 1 + 2a;^ 



(5.12) 



(5.13) 



4.t2 2.x(3 + 2a;) 4.^2(1 + 2.x) 2.^2(3 + 22;). 

The approximants p,P are not accurate for small x and tend to Too, respectively, for x — » 0+. 
For this reason, we match them with the Taylor bounds (1.29-1.30) of order m = A, which have 
the form 

t{x)<N(\,\,x)<T(x) (5J4) 

^, , 37r TT 37r 2 s 77r 4 

Tlx) := X H a;'' -| x^. 

^ ' 8 16 256 512 24576 

The final bounds for A'^ are 

1 3 

i{x) := max(t(a;),p(a;)) < N{-,-,x) < L{x) := min(T(a;), P(a;)) for x G [0,-Foo). (5.16) 

One finds numerically that: £{x) = t{x) for < a; < 1.95 and £{x) = p{x) for x > 1.96; L{x) = T{x) 
for < a; < 2.41 and L{x) = P{x) for x > 2.42. Let us introduce the relative uncertainty 
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Then: ■i?(0) = 0, '!?(!) < 0.00028; sup^g^ +00) "^i^) < 0.0074, the sup being attained at the point 
xo e (2.41, 2.42) where T{xo) = Pixo); i?(3) < 0.0051, '!?(7) < 0.0011, ^{x) ^ for a; ^ +00. o 

iii) We assume a = 2, 5 = —1/2 and apply Prop. 4. 2; in this case, to get satisfactory results it is 
convenient to increase the order of the approximation, so we choose n ~ A. In this way, we obtain 

g{x)<N{2,-^,x)<G{x) for all a; e [0, +00) , (5.18) 

where 

/ ^ 7(2,3:) 7(3, a:) 37(4,3:) 57(5,3;) 357(6,3;) 
' ■ a;2 ^ 2x^ 8x^ 16x^ 128x^ ' 

7(2, x) , 7(3, a;) , 37(4,3;) , 57(5,3:) , 4557(6, x) 

667 1087 
.(0) = ^, G(0) = ^ 

All the above incomplete Gammas are elementary: the ones of orders 2,3,4 arc read from Example 
i), the other ones are 7(5, a;) = 24 - (x" + 4a;3 + 123;^ + 243; + 24)6"^, 7(6, x) = 120 - (a;^ + 53;-* + 
20a;3 + 60a;2 + I2O3; + 120)6"^. Thus 

J_ l_ _9_ _^_/^_^ 735 399 837 645 525 \ _^ 

1 1 9 15 6825 _ / 735 2835 2499 7137 6945 6825 \ _^ 
+ ^ + ^ + ^ + (^Ym^T28^^32^^32^^ 16^^!^ • 

The relative difference 

m ■■= ^14^ (5.21) 

G{x) + g{x) 

is a decreasing function. It is: C(O) = 210/877 < 0.24, $(1) < 0.22, ^(3) < 0.15, ^(7) < 0.046 and 
^(3;) ^ for X — > +00. 

Again, for small x we can improve the previous estimates using the Taylor bounds (1.29-1.30); 
with m = 4, these are 

t{x) < N{2, -^,x)< T{x) for x e [0, +00), (5.22) 

/ X 4 16 16 2 128 3 64 4 ^ 

*(^)^=3-T5^+35^ -945^ + 2079^^ ' ^''''^ 
X 4 16 16 2 128 , 64 , 
^(^^^=3-15^+35^ -945^ + 2079^ • 
The final bounds arising from (5.18) and (5.22) are 

j{x) := max{t{x),g{x)) < N{2,-^,x) < J{x) min{T {x) , G{x)) (5.24) 

for all X G [0, +00). One finds numerically that: j{x) = t{x) for < 3: < 1.57 and j{x) = g(x) for 
X > 1.58; J{x) = T{x) for < a; < 1.54 and J{x) = G{x) for x > 1.55. The relative uncertainty 

has the following features: 77(0) = 0, 77(1) < 0.016; supj.gjQ^^) rj{x) < 0.20, the sup being attained 
at the point xq & (1.57, 1.58) where t{xo) = g{xo); r]{x) = ^{x) for x > 1.58. <> 
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Figure 4: Graphs of g{x), N{2,3/2,x) and G{x) from Example i), for x & [0, 7]. 




Figure 5: Graphs of g{x), N{2,3/2,x) and G{x) from Example i), for x e [0, 1]. 
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Figure 6: Graphs of g{x), N{2,3/2,x) and G{x) from Example i), for x € [6, 7]. 
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